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ABSTRACT 

We present matched-resolution maps of HI and CO emission in the Virgo Cluster spiral NGC 4647. 
The galaxy shows a mild kinematic disturbance in which one side of the rotation curve flattens but the 
other side continues to rise. This kinematic asymmetry is coupled with a dramatic asymmetry in the 
molecular gas distribution but not in the atomic gas. An analysis of the gas column densities and the 
interstellar pressure suggests that the H2 /HI surface density ratio on the east side of the galaxy is three 
times higher than expected from the hydrostatic pressure contributed by the mass of the stellar disk. 
We discuss the probable effects of ram pressure, gravitational interactions, and asymmetric potentials 
on the interstellar medium and suggest it is likely that a m = 1 perturbation in the gravitational 
potential could be responsible for all of the galaxy's features. Kinematic disturbances of the type seen 
here are common, but the curious thing about NGC 4647 is that the molecular distribution appears 
more disturbed than the HI distribution. Thus it is the combination of the two gas phases that 
provides such interesting insight into the galaxy's history and into models of the interstellar medium. 
Subject headings: galaxies: individual (NGC 4647) — galaxies: spiral — galaxies: ISM — galaxies: 
kinematics and dynamics — galaxies: evolution — ISM: molecules 



1. INTRODUCTION 

Galaxies that are rich in cold gas tend to have roughly equal amounts of gas in the atomic and in the molecular 
phase. But, to the best of our knowledge, star formation only occurs in the molecular phase. Therefore, the global 
balance between atomic and molecular gas is crucial to the long-term evolution of galaxies and their stellar populations. 
What determines the balance between atomic and molecular gas in a galaxy, and how does that balance evolve as the 
environment of the galaxy changes? 

On theoretical grounds we would expect that the density and temperature of the gas, the stre ngth of the dissociating 
UV field, and the metallicity all affect the relative amounts of molecular and atomic gas l)Hidaka fc Sofuel l2002t 
^Imcerecn 1993). The gas should also be in hydrostatic equilibrium with a pressure that is mostly determined by 
the stellar distribution, since stars contribute the bulk of the mass in the inner parts of galaxies where the H2 — HI 
transition occurs. 

If a galaxy is disturbed by gravitational interactio ns or by falling into an intra cluster medium, the balance between 
atomic and molecular gas may be altered. Indeed, iMiller fc OwenI l|2003j) and iMiller. Owen, fc Hill l|2003fl observe 
increased star formation activity in the members of a galaxy group or subcluster that recently merged into a larger 
cluster. In a bottom- up (cold dark matter) cosmology such disturbances should be common, and they may be important 
both to the star formation histories and the morphological evolution of galaxies. In turn, the disturbances and 
asymmetries in galaxies can inform our understanding of cosmology as they allow us to place constraints on the 
merger/interaction rate. 

The semi-analytic simulations of galaxy formation fe.g. iSomerville fc Primac3 1)19991 ): iKauffmann et alJ (|199,'^ and 

successors) are attempts to understand the observed properties of galaxies by predicting their star formation histories 
in a cosmological context. To this end the semi-analytic simulations typically assume that the star formation rate can 
be calculated from empirically motivated recipes involving the gas content and the recent merger history. However, 
these empirical recipes may not capture important phases of galaxy evolution. Thus it is desirable to have a more 
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physically motivated understanding of the interstellar medium (ISM) and star formation that can be applied to both 
s emi-empirical studies and fu ture detailed N-body/hydrodynamic simulations. 

iKrumholz fc McKed l)2005l) have recently put together a detailed model for the star formation rate in galaxies. They 
begin with the molecular cloud microphysics, in which fractal turbulence drives a small fraction of the gas to high 
enough density that it collapses to form stars. They extend the physics to galaxy- wide scales by then considering 
the molecular clouds as gravitationa lly bound entities formed in a disk that i s marginally stable against gravitational 
collapse using the Toomre criterion l)Toomre|ll964l iBinnev fc Tremainelll987j) . We note that this model still depends 
critically on knowing how much of the gas disk is in molecular form and how much is atomic. In the abs ence of 
observations, that information must be obtained from a hypothesis such as that of lBlitz fc Rosolowskvl (j2004D . which 
predicts atomic and molecular column density ratios as a function of the midplane hydrostatic pressure in a disk. 
Future work on simulations of galaxy evolution could be made more realistic by incorporation of such physically-based 
models of a multiple-phase ISM and its effects on the star formation rate. 

Our current contribution to this effort is an observational study that begins to test such physical models on a gently 
disturbed galaxy. We present an analysis of the kinematics of the neutral interstellar medium in the Virgo Cluster 
spiral NGC 4647, along with an analysis of the relationships between atomic gas, molecular gas, and star formation 
in the galaxy. We investigate whether the observed molecular /atomic balance and the disturbed kinematics can be 
understood in the context of the galaxy's environment. The form and magnitude of the kinematic disturbance in NGC 
4647 are common, so the processes affecting this galaxy might be applicable to many spirals. However, few mildly 
disturbed spiral galaxies have been mapped at matching resolutions in both atomic and molecular gas, so opportunities 
for this type of analysis are rare. 

2. MEET NGC 4647 

NGC 4647 is a gas-rich Virgo Cluster spiral whose interstellar medium has been gently disturbed and whose sur- 
roundings suggest some possible causes of the disturbance. It is 3.2° = 0.95 Mp c in pro iectiori from M87. The galaxy's 
distance has been estimated via the Tully-Fisher method to be 17. e ijGavazzi et alJlT999tl . so that it is close 
to the midplane of the Virgo Cluster. It is also notable for being in a close pair with the elliptical galaxy NGC 4649 
(M60); the two are 2.6' apart or 13 kpc in projection (Figure^). The galaxies are close in radial velocity as well, with 
NGC 4647 at 1415 ±4 km s^^ and NGC 4649 at 1114 ±17 km s"! l|de Vaucouleurs et al.. 19m.). Surface bright ness 
fluctuations give the distance to NGC 4649 to be 16.8 ±1.2 (random) ±1.2 (systematic) Mpc ()Tonrv et alJl200H) . so 
this independent distance measurement is also consistent with the two galaxies being physically close to each other. 
We therefore adopt a distance of 17 Mpc to both galaxies and entertain the possibility that they are gravitationally 
interacting. 

The location of NGC 4647 in the midst of the Virgo Cluster also suggests that the galax y may have suff ered ram 
pressure stripping by the intracluster medium. The galaxy is mildly Hl-deficient; according to lKennev fc Youn g (198^) 
the HI mass is a factor of 3 lower than what one would expect for an isolated galaxy of the same optical diameter. 
This deficiency is very typical of late- type spirals that ar e 3° in projection from M87. The molecular/atomic mass 
ratio is also typical for galaxies of the same HI deficiency fKe nnev fc Yo ung 1989). 

The stellar distribution of NGC 4647 (traced by an R-band image) and Ha emission are both modestly asymmetric, 
with an enhancement in Ha surfa ce brightness at the outer edge of the star-forming disk on the side closer to NGC 
4649 (jKoopmann fc Kennevll200^ . Here we look in greater detail at the gas distribution and kinematics to search for 
clues to the processes responsible for disturbing the ISM in NGC 4647. 

3. OBSERVATIONS AND DATA REDUCTION 

3.1. CO imaging 

^^CO 1-0 emission from the NGC 4647-4649 system was imaged with the 10-element Berkeley-Illinois-Maryland 
Association (BIMA) millimeter interferometer at Hat Creek, CA (Welch et al. 1996). The BIMA observations were 
carried out in the D configuration (projected baselines 2.3 to 11 kA) for a total of 5.5 hours on 2001 June 15, 4.25 
hours on 2003 June 8, and 2 hours on 2003 June 28. Observations were made in mosaic mode using a hexagonal 
grid of 7 pointings with center-to-center separations of 60". The 2001 observations were made with a pointing center 
of 12''43'"35.0^ +11°34'00"(J2000), and the 2003 observations were made with a pointing center of 12''43'"32.0^ 
+11°35'07". This coverage was designed to give sensitivity to CO emission throughout th e disk of NGC 4647 as well 
as in the region between the two galaxies and the central regions of NGC 4649. ISage fc Wrobell l)1989t) reported a 
single-dish CO detection of the elliptical NGC 4649, so the mosaic observations were partly designed to investigate 
those claims. FigureQlshows the field of view of these data. The velocity range covered by the CO data is 888 to 1897 
km s~^. 

Re duction of the BIMA data was carried out using standard tasks in the MIRIAD package ijSault. Teuben. fc Wright 
11995*). Absolute flux calibration was obtained from observations of Saturn and the secondary flux calibrator 3C273. 
Phase drifts as a function of time were corrected by means of a nearby calibrator observed every 30 to 40 minutes. 
Gain variations as a function of frequency were corrected by the online passband calibration system; inspection of the 
data for 3C273 indicate that residual passband variations are on the order of 10% or less in amplitude and 2° in phase 
across the entire band. 

The calibrated visibility data were weighted by the inverse square of the system temperature, then Fourier trans- 
formed. The MIRIAD task invert performs an appropriately weighted linear combination of the data from the different 
pointings. No continuum emission was evident in the line-free channels, so no attempt at continuum subtraction was 
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made. The dirty images were deconvolved using both the Clark clean algorithm and a maximum entropy algorithm, 
but the two methods gave indistinguishable results. The resulting image cube's resolutions and noise levels are given 
in Table ^ Integrated intensity images were produced by the masking method: the deconvolved image cube was 
smoothed along both spatial and velocity axes, and the smoothed cube was clipped at about 2.5cr in absolute value. 
The clipped version of the cube was used as a mask to define a three-dimensional volume in which the emission is 
integr ated ove r velocity. This ma sking method for producing moment maps is described in greater detail by Wori3 
l)2001[) and bv iReean et al.l l)2001j) . Integrating the "moment zero" image again over the spatial directions then gives 
the total CO flux of the galaxy. In addition, a velocity field is constructed by fitting a Gaussian line profile to the CO 
spectrum at each position within the galaxy. 

The total CO flux we measure from NGC 4647 is 450 Jykms"^ wit h an uncertainty of 10-15% from absolute 
calibration; it is consistent with the previous single-dish measurement. iKennev fc Youn3 l(l988^) used the FCRAO 
telescope, with a beam FWHM of 45", and they quote a total flux of 600±120 Jyk ms~^. Within the error b ars, these 
two measurements agree. It is important to note, however, that the flux quoted bv IKennev fc Yound l|198aD has been 
multiplied up by a factor of 1.5 to account for "unobserved" gas. This correction factor is an educated guess of how 
much emission was not sampled by their 5 pointings (spaced by the FWHM) along the galaxy's optical major axis. We 
can make a more accurate comparison of the fluxes by convolving our integrated intensity map with a 45" Gaussian, 
mimicking the sensitivity pattern of the FCRAO telescope. When this is done, we predict that the FCRAO telescope 
should measure 212 Jykms"^ from the central position on the galaxy. This estimate is in excellent agreement with 
the measured value of 219 ±35 Jykms"^. We conclude that the BIMA interferometer has not missed any substantial 
component of CO flux from NGC 4647. For 17 Mpc distance and a CO-to-Il2 conversion factor of 2.0x10^° molecules 
cm~^ (K km s~^)~^ the H2 mass is then (1.0 ± 0.15) x 10^ M©. The mass of helium can be included by multiplying 
up by a factor of 1.36. 

3.2. HI imaging 

HI emission from NGC 4647 was observed with the National Radio Astronomy Observatory's Very Large Array 
(VLA)^ in its D and C configurations on 2002 January 4 and 2002 December 2 and 7. We obtained 3.7 hours on source 
in the D configuration and 7.2 hours on source in the C configuration. The system was observed with one pointing 
centered in between the two galaxies. The initial D configuration observations were made with a wide bandwidth, 6.25 
MHz, in order to cover the velocities of both NGC 4647 and NGC 4649; imaging of these data revealed that all of the 
HI emission from the system was coming from NGC 4647. Subsequent C configuration observations were made with a 
smaller 3.125 MHz bandwidth (centered on the velocity of NGC 4647) and higher velocity resolution. A nearby point 
source was observed every 30 to 45 minutes as a phase calibrator. The absolute flux scale was set by observations of 
J1331-f 305 and bandpass calibration was determined from the same source. All data calibration and image formation 
was done using standard calibration tasks in the AIPS package. Initial imaging revealed which channel ranges were 
free of HI line emission. Continuum emission was subtracted directly from the raw uv-data by making first order fits 
to the line-free channels. 

The calibrated data were Fourier transformed using several different uvdata weighting schemes chosen to enhance the 
spatial resolution or the sensitivity to large-scale structures. A low resolution HI image was made covering the entire 
6.25 MHz bandwidth using the D configuration data alone and the "natural" weighting scheme (all data points are 
assigned equal weight) . A medium resolution HI image was made by combining all the C and D configuration data with 
natural weighting. The highest reso lution HI image uses the C configuration data alone and Briggs's robust weighting 
scheme l)Briggs. Schwab, fc Sram eF 1999), as implemented in the AIPS task IMAGR with a robustness parameter 0.3. 
Final images' resolutions and sensitivities arc given in Tabled Dirty images were cleaned down to a residual level of 
approximately 1.0 times the rms noise fluctuations. Integrated intensity (moment) maps and velocity fields were made 
in a similar manner as for the CO images. 

The highest sensitivity HI cube (medium resolution, combined C-l-D configurations) contains a total of 7.0 ±0.1 
Jy km s~^ of HI emission from NGC 4647. The uncertainty in this value comes from a spread among different techniques 
for making the integrated intensity image and from different sum regions. The lowest resolution HI cube gives an 
integrated flux of 6.9 Jy kms~^, consistent with the above. The highest resolution HI cube also gives a flux consistent 
with the above (7.1 ±0.3 Jykms"^). Accounting for uncertainty in the absolute Hux calibration, we derive a total flux 
o f 7.0 ±0.2 Jykms"^ of HI from NGC 4647. At 17 Mpc, this flux corresponds to (4.8 ± 0.1) x 10^ Mq of HI. 

iHelou. Hoffman fc Salpeteii H1984j) measured the HI flux of NGC 4647 with the Arecibo telescope and found 8.2 
Jykms"^ with an error of 10% to 20% (not including absolute calibration uncertainty). Warmels (1988) used the 
WSRT to measure a total HI flux of 6.6 Jykms"^ with an error of 15% to 20%. Earlier D configuration VLA 
observations of NGC 4647 are reported bv lCavatte et alJ l)1990f) . who claim an HI flux of 10.10 ±0.44 Jykms"^. It is 
not clear why the 1983 VLA dataset would give such a different flux from the current data, but we note that our present 
observations have about 5 times better sensitivity than those 1983 data when measured at the same spatial resolution. 
Thus, we believe that our new VLA map has recovered all of the HI emission from the galaxy; it is significant to note 
that even the highest resolution HI cube still detects all of the HI emission. Figure |21 shows the integrated HI and CO 
spectra of NGC 4647. 

4. RESULTS 

^ The National Radio Astronomy Observatory is operated by Associated Universities, Inc., under cooperative agreement with the National 
Science Foundation. 
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TABLE 1 

Image Parameters 





CO 


HL high res 


HI: med res 


HI: low res 


Velocity range, km 


898 - 1897 


1093 - 1738 


1093 - 1738 


626 - 1914 


Velocity resolution, km 


10.09 


10.4 


20.8 


20.8 


Beam size, " 


14.79 X 13.38 


15.63 X 13.55 


36.00 X 30.69 


45.75 x 43.75 


Beam position angle, ° 


-14.3 


-7.8 


-2.5 


-16.9 


Linear resolution, kpc 


1.2 X 1.1 


1.3 X 1.1 


3.0 X 2.5 


3.8 x 3.6 


rms noise, mjyb~^ 


100 


0.35 


0.18 


0.43 


N(HI) sensitivity, cm~^ 


2.8x10^0 


5.7x10^^ 


l.lxlO^^ 


1.5x10^^ 


N(H2) sensitivity, cm~^ 









Note. — All velocities in this paper are in a heliocentric system using the optical definition. 
The distance to NGC 4647 is taken to be 17 Mpc. Column density sensitivities are 3a in one 
channel; the H2 limit assumes a CO-to-H2 conversion factor of 2.0x10^" molecules cm~^ per 
K km s^^ of brightness temperature in ^^CO 1-0. 



4.1. HI distribution 

Figure O shows the HI column density in NGC 4647 along with images of the CO intensity, an R-band image, and 
continuum-subtracted Ha, all at the same linear scale. We have also computed the radial distribution of HI in NGC 
4647 from azimuthal averages over elliptical annuli in the approaching and receding halves of the galaxy (Figure 0]). 
"Approaching" and "receding" are defined with respect to the adopted center of the galaxy and the kinematic major 
axis that are described in Section below. 

The HI column density in NGC 4647 is relatively flat at 7.5 Mq pc~^ in the interior of the galaxy. (Gas surface 
densities quoted in Mq pc~^ in this paper are corrected to face-on and include a factor of 1.36 for helium; column 
densities in cm~^ are as observed.) One notable feature of the HI distribution is an arm or segment of a ring with 
enhanced column densities at radii of 30-45" — 2.5 to 3.7 kpc on the east side of the galaxy. HI column densities 
are 10% to 20% higher in the arm than in the central plateau; the highest HI column densities in the galaxy are 
here. While the arm is strongest on the east and north sides of the galaxy, it is weak or nonexistent on the western 
side. Beyond the arm, on the north and east sides of the galaxy, the HI column density declines steeply; it falls below 
1.25x10^° cm~^ at a radius of 67" (5.5 kpc), just a little over one FWHM beyond the ridge of the arm. 

The second notable feature of the HI distribution in NGC 4647 is a very shallow dropoff in column density toward 
the southwest. There is no well defined ridge or arm on this side, and whereas the HI extends to 67" on the northeast 
side of the galaxy it peters out to 90" (7.4 kpc) on the southwest side. Both radii are measured at the same 1.25x10^° 
cm~^ in the data cube with 15" resolution. At low resolution (45"), the HI column density drops below the sensitivity 
limit of ~ 2 x 10^^ cm~^ at 130" (10.7 kpc) radius on the southwest side; however, there is relatively little HI at radii 
beyond 90" since the low resolution map does not contain any more flux than the hi gh resolution map (Section rO|l . 
The larger HI extent on one side is reminiscent of the asymmetric HI in MlOl (Baldwi n. Lvnden-Bell. fc Sancisll OSO). 
Curiously, though, the southwest side of NGC 4647 has the shallowest HI gradient but this is also the side with the 
steepest decline in CO intensity (Sect ion 14.211. 

In the earlier map of NGC 4647 by iCavatte et alJ (|1990) , the HI appeared to be slightly extended in the direction 
of NGC 4649. Those authors suggested that this extension could be caused by gas transfer from NGC 4647 to NGC 
4649, possibly fueling the nuclear activity in NGC 4649. Our present images show no sign of any HI extension towards 
NGC 4649, so we conclude that we are not presently observing gas transfer from the spiral to the elliptical. 

4.2. CO distribution 

The major features of the molecular distribution in NGC 4647 are a large concentration near the center of the galaxy, 
a plateau on the east and a steep radial dropoff to the southwest (Figures |2| and [SJ . 

The central molecular clump is poorly resolved by these 15" data, but thirty percent (4.1 x 10^ Mq, including helium) 
of the molecular mass of the galaxy is found within 15" of the central peak. At this radius the average circular velocity 
of the two sides of the galaxy is 77 km ('Section l4.3|) . so the implied dynamical mass is 1.7x10^ Mq or only four 
times the mass of the molecular concentration. The central molecular clump thus makes a non-negligible contribution 
to the dynamics in the interior of the galaxy. The molecular peak is also 3" (250 pc) west of the i?-band nucleus, and 
this offset may have important consequences for the gas kinematics in the nuclear regions. 

The eastern half of the galaxy shows an extended molecular disk that is coincident with the ring or arm segment 
seen in HI. While the atomic gas shows a 20% enhancement in column density from one side to the other at this 40" 
radius, the molecular gas shows a much stronger enhancement. Molecular surface densities (corrected to face-on and 
including helium) are 30 Mq pc~^ in the arm/ring and are a factor of two to three smaller on the west side (Figure 
ISJ. Thus the CO emission is strongly asymmetric and extends twice as far from the galaxy center on the eastern side 
as on the western side. 

Because of the asymmetry in the CO distribution, the western side of the galaxy shows a gradual decrease in atomic 
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gas surface density but a very sharp cutoff in molecular surface density. The molecular distribution of the galaxy then 
appears to have an edge at a position angle of roughly —28°. This cutoff also makes itself visible in the spectra (Figure 
13), where it is plain that the CO does not extend far enough on the low velocity side of the galaxy to sample the flat 
part of the rotation curve. The HI spectrum shows both "horns" but the CO spectrum shows only one horn. 

Star formation activity in NGC 4647 roughly follows the CO surface brightness. Ha intensities are highest in the 
nucleus of the galaxy, on the central CO clump, and are lowest on the western side where there is little molecular 
gas (Figure 13). A few small HII regions on the western side suggest that there is probably a modest amount of 
molecular gas there but when averaged over the 15" beam size it is below our sensitivity limits. We also estimate the 
star formation efRciency (SFE) in the molecular gas by convolving the continuum-subtracted Ha image to the same 
resolution as the CO image and dividing the two. By this measure the SFE is not significantly different on the two 
sides of the galaxy (at least over the regions where molecular gas is detected) . The SFE in the galaxy nucleus is also 
comparable to that in the disk, but the southeast quarter of the galaxy has a SFE a factor of two to three higher than 
the northeast quarter. Local molecular peaks in the arm of the galaxy tend to occur in between HII regions, when 
measured on the scale of our CO resolution (1.2 kpc). It would undoubtedly be useful to investigate possible effects 
of extinction on the Ha surface brightness through FIR or radio continuum images. 

4.3. Gas Kinematics 
4.3.1. Asymmetry 

The atomic gas in NGC 4647 shows a striking asymmetry in its kinematics (Figure|Sl). In the western (approaching) 
half, the side more distant from NGC 4649, the rotation curve rises quickly and flattens, giving a typical spider-diagram 
pattern with isovelocity contours that bend toward the major axis. However, on the eastern side of the galaxy the 
velocity contours are much straighter and more suggestive of a linearly rising rotation curve. This pattern is also 
evident in the individual channel maps (Figure [71), which show the bent butterfly- wing structures on the west side of 
the galaxy but linear features in the east. Similarly, the major axis position- velocity diagram (Figure [HJ shows that 
the west side rotation curve (RA offsets < 0) flattens whereas the east side kinks at a radius of 10" and continues to 
rise. 

4.3.2. Kinematic parameters 

Analysis of the HI velocity field allows us to estimate the important kinematic parameters of NGC 4647 and to 
quantify the kinematic asymmetries discussed above. We carried out a set o f tilted ring fits with the rot cur task in the 
GIPSY software package, according to the method described bv iBegemanI (^^3) and summarized bv iSwaters et alJ 
llT999.) . 

We first determine the position of the optical nucleus in NGC 4647 and compare that to the kinematic center of the 
galaxy. Ecliptic coordinates were transferred from a POSSII Digitized Sky Survey plate to the nearly-Mould-R image 
of KooDmann. Kcnnev. & Young (2001) using seven stars distributed uniformly around NGC 4647 and NGC 4649. 
The position of the galaxy nucleus in this image is 12'' 43™ 32.58", +11° 34' 56.8" (epoch J2000) with an uncertainty 
of half a pixel (0.4"). A kinematic center is determined by first fitting all annular rings in the galaxy with all of their 
parameters (center, position angle, inclination, velocity) left free, and the average of all of the fitted centers gives a 
kinematic center position only 2.5" away from the optical position. Additionally, a very simple fit of an inclined disk 
with a rotation curve V{r) oc 1 - e''/''« (as encoded in the AIPS task GAL) gives a kinematic center 4" away from the 
optical nucleus. As the resolution of the radio data is 15", we consider the kinematic centers to be consistent with the 
position of the optical nucleus and we adopt the latter for the dynamical center of the galaxy. 

Subsequent tilted ring fits with fixed centers show significant trends in the systemic velocity from 1407 km s^^ at 
radii less than 20" to 1425 km s~^ at 60-80". We attribute this trend to the asymmetry in the velocity field - more 
specifically, to the fact that the approaching side rotation curve flattens whereas the receding side rotation curve keeps 
rising. A value of 1410 km s^^ is adopted for the systemic velocity since it minimizes the differences between the two 
sides at radii < 25". 

The kinematic position angle of the galaxy is significantly better constrained than the systemic velocity. Tilted ring 
fits for position angle (with kinematic center fixed as described above) give position angles 98°- 99°, and no trend with 
radius. We also fit the model exponential velocity field to receding and approaching sides independently and find that 
both sides of the galaxy give angles consistent with the tilted ring fits. We take the position angle of the kinematic 
major axis to be 98.5°± 0.5°, where the uncertainty indicates th e dispersion I j etwee n various fits. Note, however, that 
this position angle is significantly different from that adopted bv lRubin et alJ l)1999|) for longslit Ha kine matics, which 
was 124°. It is now clear that the reason NGC 4647 appears to show minor axis rotation in the work of 'Rubin et alJ 
(^999) is that their adopted kinematic minor axis is 30° offset from the true orientation- we observe no minor axis 
rotation when the proper kinematic position angles are used. 

The kinematic major axis obtained from the HI velocity field is in good agreement with the near-IR morphological 
major axis. We obtained the Two Micron All-Sky Survey (2MASS) K-b and image of the NGC 4649-4647 system, 
masked out NGC 4647, and used the Multi-Gaussian Expansion routines of lCappellari ) f2002) to model NGC 4649 and 
remove it. In subsequent MGE fits of NGC 4647 itself and isophote fits using the ELLIPSE routine in STSDAS, we 
find the isophotal major axis position angles between 97.6° and 100.7°. 

The receding half of the galaxy, with its almost linear isovelocity contours, gives very poor constraints on the 
inclination of the disk. Therefore the best information on the inclination comes from fits to the approaching half alone. 
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with fixed center and position angle. Tilted ring fits are consistent with the fit of the simple exponential velocity field 
and b oth give inclinations 39.3° ± 1.8°. The value is also consistent with the inclination adopted bv iR^ibin et alJ 
l)1999|) as determined from optical isophotes. Figure IHl shows the final rotation curve derived from the tilted ring fits, 
where the receding and approaching sides were fit independently (weighted by the square of the cosine of the angle to 
the kinematic major axis). 

The magnitude and spatial distribution of the kinematic asymmetry in NGC 4647 can be estimated by taking the 
rotation curve derived from the approaching (flatter) side alone, constructing a model velocity field, and subtracting 
from the observed velocity field. The residual velocities are on the order of 20 km s~^ at r ~ 60" on the east side of the 
galaxy. If these residual velocities arise in tangential motions in the plane of the disk, then correcting for inclination 
gives a 30 km s^^ velocity difference from one side of the galaxy to the other (compared to a circular speed of 115 
km s~^). 

4.3.3. Atomic vs. Molecular and Ionized Gas Kinematics 

Molecular and ionized gas kinematics in NGC 4647 agree well with the atomic gas. For example, aside from the 
fact that the CO is more centrally concentrated than HI and does not extend as far as HI on the western side, the 
channel maps in Figure dare nearly identical in CO and HI. The position- velocity diagrams in Figures 151 and [TUI also 
show good agreement between HI, CO, and Ha kinematics. 

There is one notable difference between molecular and atomic gas kinematics: Figure |H1 shows that at radii between 
10" and 30" on both sides of the nucleus the HI line profiles are asymmetric. The line profiles have tails extending 
to smaller rotation speeds, which is symptomatic of gas lagging the rotation curve by as much as 40 km s^^ (60 
km s~^. after correcting for the inclinatio n of the disk). This HI "beard" has also been observed in other spirals 
ijFraternali et aljr2001t iBarbieri et al.ll2005j) and is interpreted as gas above and below the plane of the disk. We note 
that the molecular gas does not show these asymmetric line profiles, so the gas that is inferred to be up out of the 
plane is primarily atomic. 

4.3.4. Dispersions, Disk Instabilities, and Star Formation 

Gaussian fits to individual line profiles are used to estimate the velocity dispersions of the HI and CO. The con- 
tribution to the linewidth from beam smearing is estimated at each point as the local dispersion in the velocity field 
(weighted by a Gaussian of the same size as the spatial resolution). Broadening by beam smearing and finite channel 
width are then subtracted from the line width in quadrature. Aside from the nuclear regions where beam smearing 
dominates, the corrected HI dispersions in the disk are 12± 3 km s~^; the latter quantity indicates the width of the 
distribution rather than an uncertainty in measurement. The molecular gas has somewhat narrower lines, consistent 
with dispersions 7 ±3 km s~^. There is no measurable difference in the velocity dispersions on the two sides of the 
galaxy 

The star formation rate in the disk may depend not only on the gas surface d ensity, but also on the degree to which 
the gas disk is gravitationally unstable to the formation of m olecular clouds iKrumholz & McKee 200^1) . The disk 
instability is normally parametrized by Toomre's Q parameter l)Toomrelll964t iBhm^Tfc'^emaindllQS'ijr ^ 



which is based on a comparison of the local sound speed a and the epicyclic frequency k to the gas surface density 
Sgas- Large values of the Q parameter indicate stability against the formation of self-gravitating entities. 

Given the differences in the properties of the two sides of NGC 4647's disk, the galaxy makes an interesting test 
case for the study of star formation. The gas velocity dispersion is equal on the two sides of the galaxy; however, the 
epicyclic frequency may not be. Under the assumptions that the HI and CO velocity fields still trace the local circular 
velocity and that the inclination of the disk is the same on both sides of the galaxy, we may use the results in Figure 
Elto estimate the epicyclic frequency. The side with the rising rotation curve has k = 2.1 x 10~^^ s~^ at a radius of 4 
kpc, whereas the flatter side of the rotation curve has k = 1.2 x 10~^^ s~^ at the same radius. These initial estimates 
hint that the higher gas surface densities on the east side of the galaxy may be offset by a higher epicyclic frequency 
on that side, producing roughly equal values of Toomre's Q parameter on both sides. 

4.4. Ehi, 5]h2, and Pressure 

The asymmetry in the CO intensity (Section 14. 2|) could conceivably reflect either the underlying distribution of 
molecular gas or an asymmetry in the CO luminosity per unit H2 mass. For example, a factor of two to three variation 
in the CO-to-H2 conversion factor from one side of the galaxy to the other could reproduce the observations while 
keeping the total gas surface density equal on both sides. This conversion factor is expected to depend on several 
properties of the molecular clouds, including the gas density, temperature, and the interstellar UV field to which the 
clouds are exposed. Here we simply comment that the UV field is probably stronger in the eastern half of the galaxy 
than in the western half, since the Ha luminosity is high er in the east (Fig ure Et. A stronger UV field could either 
lead to higher temperatures and increase d CO luminos ity ijWeiss et alJl200l]) or to preferential destruction of CO and 
lower CO luminosity (Israel 1997; Leaue ux et all 119941 ) . and distinguishing between these alternatives is beyond the 
scope of the present paper. 

We also note that if the conversion factor varies but the total gas surface density is constant, then the star formation 
efficiency must be higher on the eastern side of the galaxy in order to produce the higher Ha intensities there (Section 
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I4.2|l . A higher star formation efficiency on the eastern side might be in conflict with the Toomre-style stabihty 
analysis of section IT. 8. 41 which showed that the epicyclic frequency is expected to be higher in the east. Evidently the 
assumption of a variable conversion factor raises at least as many questions as it answers. Thus we proceed under the 
assumption that the conversion factor is constant and the galaxy really does have factors of two to three higher total 
ga s surface de n sity on the eastern side than on the western side. 

iRhtz Rosnlowskvl lf2004 ) propose that the H2/H I column density ratio in a disk galaxy is determined by the 
hydrostatic pressure. This kind of a model could have important applications in simulations of galaxy evolution, 
since it provides a way to calculate the molecular surface density for star formation purposes. The asymmetric gas 
distribution in NGC 4647 provides an unusual opportunity to test this model. From a different perspective, if the 
model is correct then the observed column densities in NGC 4647 offer a way to calculate the interstellar pressure on 
the two sides of the galaxy and to find out which (if cither) of the two sides is consistent with the properties of other 
spirals. 

We expl ore these questions by constructing a map of the hydrostatic pressure in the galaxy using the method of 
iBlitz fc Ro solowskv ( 2006) . Their formulation models disks of galaxies as two gravitationally-coupled fluids with the 
stellar scale height and midplane volume density much larger than those of the gas at any given position in the disk. 
For an infinite disk, the implied midplane pressure in the gas component is 

^ = 272 cm-3 K { \ { T { { " ^ • (1) 

fc V^^0 PC-V VA^0 PC-V Vkms-VVPC/ 

where Sg and E^, arc the inclination corrected surface densities of neutral gas and stars respectively. The gas velocity 
dispersion is given by Vg and h^, is the stellar scale height. We convert the integrated intensity of the H I and CO 
maps into mass surface densities using the standard conversion factor of 1.82 x 10 cm^ -V(K km s"^) for the 21 -cm 
emission and a CO-to-H2 conversion factor of 2 x 10^'^ cm^^/(K km s^^). The column densities are then converted 
to mass surface density assuming a mean particle mass of 1.36 toh per hydrogen nucleus and an inclination of 39.3° 
f ti4.3.2|l . We convert the masked j^T-b and image of the galaxy to a stellar surface density assuming a mass-to-light 
ratio of Mk/Lk = 0.5 Mq/Lq l|Bell fc de JonEil200lj) . After convolving the H2 and stellar surface density maps to 
match the resolution of the H I map, we aligned the H I and stellar surface density maps to the H2 map. Using the 
resulting maps of and Eg ~ Ehi + Sh2, the external pressure can be c alculated directly fro i n Equ ation [H provided 
the gas velocity dispersion and stellar scale height are known. Following "Blitz fc Rosolo wskv! l|2006|) . we assume that 
a gas velocity dispersion of 8 km s~^ appropriate for most galaxie s and determine a stellar scale-height (420 pc) from 
a fit to the stellar scale-length (3.2 kpc) using the relationship of fKregel et al.l 1|2002D . 

In Figure 1111 we plot the molecular gas fraction as a function of the midplane hydrostatic pressure Pext iu the 
two halves of the galaxy (separated along the kinematic minor axis). For the western half of the galaxy, the surface 
density ratio scales with the hydrostatic pressure as it does in other galaxies (Eh2/Shi = {Pext/ PoY'^"^^'^)- The scale 
con stant Pn (= 1-5 x 10"* K cm~'^) is consistent with the other galaxies in the Virgo Cluster that were included in 
the IBlitz fc R.o solowskv (2006) sample. In the eastern half of the galaxy, the surface density ratio for low apparent 
pressures (large radii) is a factor of 3.3 higher than in the western half. In other words, the eastern half is proportionally 
more molecule-rich than the western half for a given apparent pressure in the ISM. Provided that the pressure in the 
ISM does indeed regulate the molecular gas fraction, this result implies that the eastern half of the galaxy must be 
subject to larger pressures in the ISM than are inferred from the mass components of NGC 4647 alone. The pressure 
deficit is P/k — 2.4 x 10** K cm"'^, and this deficit is twice as large as the midplane hydrostatic pressure in the outer 
disk of the western side of the galaxy. 

5. DISCUSSION 

In a ACDM universe galaxy-galaxy mergers and encounters are recognized as major processes that have determined 
the spectrum of galaxy properties. From an observational perspective, galaxy asymmetries can serve as indicators 
of the merger/interaction rate. One may measure the incidence of asymmetries and then make the leap to the 
merger/interaction rate via several important assumptions such as the timescale over which asymmetries will remain 
detectable. That timescale depends crucially on an understanding of the processes responsible for the asymmetries 
and how they affect the gas, stars, and star formation activity. 

In this context our observations of the mildly disturbed galaxy NGC 4647 raise several questions. NGC 4647 
shows two kinds of lopsidedness, in its kinematics and in its gas distribution. More specifically, the gas distribution 
is asymmetric in two ways: the total gas surface density is significantly higher on the east side than on the west, 
and the gas disk is also more highly molecular on the east side than on the west. How common are these types 
of asymmetries, and what causes them? We address the problems in two stages. In section 15.11 a comparison to 
other mildly lopsided galaxies suggests processes that may create the kinematic asymmetry and enhanced gas surface 
densities on the east side of the galaxy (or the depressed surface densities on the west side). Section l5?^ discusses 
local pressure enhancements that may drive the ISM into a more highly molecular state on the east side of the galaxy. 
In section 1^31 we list some remaining questions and propose tests that would give greater insight into NGC 4647 in 
particular and spirals in general. 



5.1. Lopsided Galaxies 
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It h as bee n known for some time that many spirals are asymmetric or fopsided. For example, iRix fc Zaritskvl 
lll995D and iZaritskv & Rixl (119971^ fi nd that some 30% of field spirals have lopsided ste llar distributions. 
iBaldwin. Lvnde n-Bell. fc Sancisi* lll980l) show asymmetries in the HI distributions of spirals, and iRichter fc Sancisil 
l)1994|) find that 50% of field spirals show asymmetric global HI profiles that are probably attributable to lopsided gas 
distributions. Rubi n et alJ ((1999) studied Ha rotation curves for Virgo spirals, and they find that half show significant 
kinematic disturbances. In many of these, the disturbance takes the form of an asymmetry in which the velocity field is 
different on the two sides of the galaxy. From studies of HI position- velocity diagrams and velocity fields lSwaters et all 
l)1999|) estimate that the incidence of kinematic lopsidedness is 15%-30% among non-interacting field spirals. 

The relatively high incidence of asymmetries (even among field spirals) is significant, since it implies that the re- 
sponsible mechanisms must either be frequent or long-lasting. Typical dynamical timescales in spiral disks are only on 
the order of 10^ yr. In this context, four general processes have been proposed as sources of kinematic and/or morpho- 
logical lopsidedness: ram pressure in a cluster environment, gravitational interactions between galaxies, asymmetric 
accretion of gas from the cosmic web or from satellites, and asymmetric or off-center gravitational potentials. 

NGC 4647 is undoubtedly experiencing some ram pressure in its passage through the Virgo intracluster medium. 
Section 15.21 makes quantitative estimates of the ram pressure. At a purely qualitative level, the sharp HI edge on 
the northeast side and the shallow dropoff in the HI column density on the southwest fSection l4.1|) are suggestive of 
gentle stripping. The kinem atic asymmetry in NG C 4647 may also be evidence of the non-circular velocities induced 
by significant ram pressure ijHidaka fc Sofuell2002(l . Finally, the form of th e CO asymmetry in NGC 4647 is very like 
that in NGC 4419, which was attributed to ram pressure bv lKennev eiTaLl Cl990). 

However, it is not at all obvious that ram pressure by itself is a convincing explanation for the ensemble of properties 
displayed in NGC 4647. The east side of the galaxy appears strongly disturbed, with its enhanced CO emission and 
rising rotation curve. At 40" radius the total gas density is twice as high on the east side as on the west side. In 
contrast, the west side appears very gently (if at all) disturbed. Its rotation curve is well behaved out to at least 70". 
Unless the galaxy is proceeding exactly edge-on through the intracluster medium, so that the west side is shielded 
by the east side, we might expect that the effects of ram pressure would be more noticeable on the low density side 
of the galaxy than on the high density side. In this respect the ram-pressure stripped galaxy NGC 4419 is actually 
quite different from NGC 464 7; both have asymmetric molecular dis tributions, but in NGC 4419 the atomic gas is 
very severely stripped as well l)Kennev et aLllloO^ iChung et al.ll2006|) . 

Asymmetric gravit ational potentials can also produce kinem atic features like those observed in NGC 4647 and other 
spirals. For examole. lSchoenmakers. Franx. fc de Zeeuwl l)1997r add a perturbation of the form cos(m6') to the potential 
of an axisymmetric galaxy and calculate the effects on its velocity field. Based on that work Swaters ct al. ( 1999) have 
shown that a to = 1 perturbation in the potential, where the axis of the perturbation is aligned with the kinematic 
major axis of the galaxy, can produce a velocity field that fiattens on one side but continues to rise on the other. A 
5%-10% perturbation in the potential can induce 10%-20% amplitude differenc es in the rotation cur ves of the two 
sides of the galaxy. This kind of a lopsided potential is the explanation favored bv lSwaters et a l.* fl999') to explain the 
HI velocity fields of DDO 9 and NGC 4395, which are similar in character to NGC 4647. iNoordcrmccr ct al. (200l|) 
have also shown that the lopsided kinematics of NGC 4395 (and, by analogy, NGC 4647 ) could be explained by setting 
the galaxy disk off-center with respect to the dark matter halo. iBattaglia et al.l l|2006j) also propose that the halo of 
NGC 5055 is offset from the center of its disk. 

A lopsided gravitational potential would affect the stellar distribution in NGC 4647 in addition to the gas distribution. 
There are subtle hints of a to = 1 asymmetry in the stellar distribution of the galaxy; for example, in the broadband 
red image of Figure 13 the outermost contours are compressed on the southern side of the galaxy and extended on the 
northwest. The galaxy nucleus is therefore not quite centered in the middle of the outer contours. NGC 4649 has been 
modeled and subtracted, but more careful analysis should be done to quantify any degree of asymmetry in the stellar 
distribution of NGC 4647. 

NG C 4647, of course, is also an obvious candidate for a gravitational interaction with NGC 4649. iBournaud et alJ 
l)2005l) have modelled distant, high velocity encounters of similar-mass galaxies, and they find that for impact param- 
eters 130-450 kpc and velocities 160-450 km s~^, the retrograde in-plane encounters can make m = 1 asymmetries 
in which the ratio of the perturbed to the unperturbed gravitat ional force is greater than 0.1 for several Gyr. These 
amplitudes are similar to those required bv lSwaters et al.l l|1999l ^. 

A statistical analysis of the distribution of galaxy lopsidednesses leads IBournaud et al.l ()2005t) to propose that 
asymmetric accretion of gas from the cosmic web must be a significant contributor to lopsidedness in general. The 
mechanism simply requires that the incoming gas have a nonzero impact parameter. They find that if accretion rates 
are fairly high — several Mq yr~^, or enough to double the galaxy masses over a Hubble time — a galaxy may develop 
an to = 1 asymmetry with an amplitude of 10% to 20%. The asymmetries in the gas and stellar distributions can 
persist even for several Gyr after accretion stops. As NGC 4647 is embedded in the Virgo cluster's hot ICM it is 
unclear how this process might apply to NGC 4647, what the phase of the accreted gas would be, or whether the effect 
on galaxy kinematics would match the observations. 

We suspect that the process that is responsible for the kinematic asymmetry also enhances the gas densities on 
the eastern side of the galaxy. But the opposite perspective is also interesting: the pile-up in gas surface density 
may contribute to the kinematic asymmetry. We estimate the dynamical mass of the galaxy taking a rotation speed 
of 115 km at a radius of 70" on the approaching (flatter, probably less disturbed) side of the rotation curve; 
these values give 1.8x10^" Mq within 5.8 kpc. Some 7.7% of this dynamical mass is attributable to molecular gas 
plus helium. Since most of the molecular gas is found on the east side of the galaxy, that gas distribution may well 
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contribute to a few percent perturbation in the gravitational potential. The gas density asymmetry is also a m = 1 
perturbation whose axis is roughly aligned with the kinematic major axis, so it has the correct form to explain the 
kinematic asymmetry according to the models of lSwaters et al.l l)1999|) . Therefore, even if the original cause produced 
both the kinematic and the gas density asymmetries together, the resulting density asymmetry may help to maintain 
the kinematic asymmetry over timescales longer than a dynamical time. 

5.2. Producing the Pressure Asymmetry 

If pressure regulates the chemical state of the ISM, as suggested bv lBlitz fc Rosolowskvl l)2006j) . then the asymmetry 
in molecular gas fraction is produced by an asymmetry in ISM pressure in the disk. Pressure, by itself, cannot explain 
the enhanced gas column density on the east side of the galaxy. However, once the gas is piled up on one side, pressure 
could explain the amount of the ISM that is molecular. The analysis of section ^31 finds that the observed asymmetry 
in the molecular gas fraction would imply an excess pressure oi P/k ~ 2.4 x 10'' K cm^'^ beyond that required for 
the hydrostatic support of the disk of NGC 4647. Thus, we consider again the effects of the intracluster medium, the 
proximity of NGC 4649, and asymmetric gravitational potentials as possible sources of extra pressure on the east side 
of NGC 4647. 

At the ^ 1 Mpc separation between M87 and NG C 4647, the intercluster medium (ICM) of Virgo has a density 
of 71 10^^ cm^'^ and a temperature of ~ 2 keV ijKennev &: Youn3ll989j) . The isotropic pressure in the ICM is 
significant: Picm /k ~ 2 x 10"^ cm~^K, on the order of 10% of the midplane hydrostatic pressure in the outer parts of 
NGC 4647's disk. If NGC 4647 is moving with respect to the ICM, it suffers additional ram pressure 

Pram = = 1.3 X 10^ cm-3 K f ^^^/ . (2) 

z \800 km s / 

The equation has been normalized to the one-dimensional velocity dispersion of the Virgo Cluster ijHuchral 1985(1 . NGC 
4647 has a line-of-sight motion relative to M87 of only 100 km s^^, but if its tangential velocity carries it eastward on 
the sky at 1000 km s~^ or so then the ram pressure could be strong enough to effect the inferred pressure difference. 
It remains to be seen whether the extra pressure could be applied delicately enough to avoid major disruptions in the 
kinematics of the downstream side of the galaxy. 

Provided that NGC 4647 is physically close to NGC 4649, the hot gas in NGC 4649 can increase the pressure on 
the eastern side of NGC 4647 through two mechanisms. The first is ram pressure between NGC 4647 and the halo of 
NGC 4649. Indeed, the diffuse X-ray emission fro m NGC 4649 extends out to the projected distance of NGC 4647 
with a suggestion of a bow shock between the two ijRandall et a l. 2004). If the galaxies are approaching each other, a 
lower limit on the relative velocities of the gas systems is 1^647 + K ot cosi — ^4649 =400 km s~'. An enhancement in 
the plasma density of 0.002 particles cm~'^ on the eastern side of NGC 4647 would produce thejiecessary pressure to 
accou nt for the observed asymmetry. Using the density profile of the X-ray halo of NGC 4649 (Brighcnti & Mathc^^ 
Il997() . we estimate that this particle density occurs at roughly 12 kpc from the center of NGC 4649. The ISM of NGC 
4647 may also be compressed by the isotropic pressure of the hot halo gas; again, the thermal press ure of NGC 4649 is 
2.4 X 10** cm"^ at a distance of ^ 12 kpc from the center of the galaxy |Brighcnti & Math ewjll997f) . The two galaxies 
would have to be no farther apart than their current projected separation for these mechanisms to drive the molecular 
asymmetry. 

The difficulty with appealing to a close approach of NGC 4647 and NGC 4649 to explain the kinematic and pressure 
asymmetries is, again, it seems unlikely that the disk of NGC 4647 could survive the tidal forces and remain unscathed. 
We take the m ass of NGC 4649 to be at least 10^^ Mq based on the total V mag nitude and a F-ba i id ma ss-to-light 
ratio of 9.5-16 l)De Bruvne et alJl200lD . Strictly speaking, the dynamical analysis of lDe Briivne et all l)200lD measures 
the mass-to- light ratio within the effective radius so that this dynamical mass of 10^^ Mq is probably an underestimate. 
For galaxy separations as small as 30 kpc, the tidal acceleration at the edge of NGC 4647's HI disk is two-thirds as 
large as the spiral's own centripetal acceleration. We therefore consider it extremely unlikely that NGC 4647 and 
NGC 4649 are presently receding from a past encounter at a few tens of kpc. It may be possible that they are still 
approaching each other. 

Alternatively, if something (NGC 4649, or an off-center dark matter halo) generates a mild m = 1 mode in the 
potential of NGC 4647 as described in section 15.11 the gas must speed up and slow down as it executes its nearly 
circular orbit. The ram pressure of gas flowing down into the trough in the potential may be sufhcient to account for the 
expected asymmetry in pressure. This possibility is especially appealing since it naturally connects the asymmetries 
in the velocities, the gas surface density, and the molecular fraction. The ram pressure of such a flow is 

Pram = ^M^^ = 3 X 10^ cm-^ K ( ^ (3) 

2 V cm / \ 10 km s / 

where Vfiow is the relative velocity at which gas approaches the potential minimum. For flow velocities comparable 
to the observed velocity differences on the two sides of the galaxy, the combination of the ram pressure with the ICM 
and the pressure of the gas accumulating in the m — 1 potential could explain the observed gas distributions. 

Finally, since the implied hydrostatic pressure scales linearly with the gas velocity dispersion, we note that the 
molecular asymmetry could be explained without recourse to an additional source of pressure if the gas velocity 
dispersion is higher on the east side of the galaxy than on the west. However, to account for the molecular asymmetry 
entirely, the velocity dispers ion wo uld have to be > 20 km s~' on the east side. Such high dispersions are ruled out 
by our observations fsection [4.3.4|l . 
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5.3. Remaining Questions 

How long the asymmetries in the kinematics and the gas distribution of NGC 4647 wiU last depends critically on 
their causes. If the disturbance is temporary, such as a fly-by of NGC 4649, then the spiral can be expected to settle 
back to a more symmetric state eventually. Alternatively, an offset dark matter halo might produce a long-lasting 
asymmetric state. The discussions above do not provide a definitive answer to the source of the asymmetries in NGC 
4647, but a variety of observational and simulation-based tests can help confirm or rule out various processes. 

Part of the difliculty with understanding NGC 4647 is that the separation between NGC 4647 and NGC 4649 is 
still unknown (as is the true relative velocity of the two). However, if the asymmetries of NGC 4647 are due to a 
galaxy-galaxy interaction then we would expect the stellar kinematics of NGC 4649 to refiect that interaction. N-body 
simulations could be used to predict the character and the magnitude of those interaction signatures. 

We have also discussed whether hydrodynamic processes such as ram pressure stripping or off-center gas accretion 
could be responsible for the kinematic and gas distribution features of NGC 4647, but we do not know the history of 
NGC 4647's interaction with the ICM. The importance of hydrodynamic processes could be assessed by comparing 
the gas kinematics to the kinematics of the old stellar populations. If the stellar kinematics are symmetric, then the 
asymmetries in the gas can probably be attributed to temporary hydrodynamic effects. In addition, hydrodynamic 
simulations could be used to test whether the ensem ble of pro perties displayed by NGC 4647 implies the presence of 
both ram pressure and gravitational disturbances (as iVoUmerl (2003) has suggested fo r NGC 4654). 

A comparison of the m = 1 asymmetry in the gas surface density to the models of iSwaters et al.l l|1999D suggested 
that the lopsided gas distribution might contribute to the observed kinematic asymmetry in the gas. That idea could 
be tested with appropriate hydrodynamic simulations. In addition, simulations that can follow heating, cooling, and 
formation and dissociation of molecules could investigate whether the conversion from HI into H2 follows naturally 
from the differences in the gas velocity from one side to the other. 

Ha images of NGC 4647 are noticeably asymmetric (Figure IS)), and we have shown that the Ha surface brightness 
roughly traces the CO surface brightness. If the H2/CO conversion factor is roughly constant throughout the galaxy, 
the star formation efficiency is constant as well. Yet the epicyclic frequency, an important component of a local 
gravitational stability analysis, may be significantly different on the two sides of the galaxy. More detailed analyses of 
the stability of this asym metric disk could give valu able insight into large-scale star formation processes. 

The pressure model of iBlitz fc R.osolowskvl l)2004j) is consistent with our observations of NGC 4647, to the extent 
that there are several viable candidate sources for the implied extra pressure on the east side of the galaxy. Further 
progress in testing the model could be stimulated by a better knowledge of how many spirals show relatively symmetric 
atomic gas but strongly asymmetric molecular gas. 

We note that NGC 4647 is actually classified as kinematically regular (non-disturbed) bv lRubin et alJ l|1999D . Its 
kinematic asymmetries are indeed milder than those of many Virgo spirals, as one can see by perusing the rotation 
curves presented by those authors. But the asymmetries of NGC 4647 are more obvious in two-dimensional velocity 
field data than in one-dimensional longslit data. The properties of NGC 4 647 therefore sugge st that the true incidence 
of disturbances among Virgo spirals is higher than the 50% measured bv lRubin et alJ l)1999t) . 

Furthermore, NGC 4647 would probably not have been identified as asymmetric on the basis of its global HI profile 
alone (Figure I^J. The galaxy therefore suggests that the incidence of asymmetric gas distributions may be higher than 
the 50% rate inferred by Richter & Sancisi (1994). It is only the combination of the HI and the CO data that revealed 
the true peculiarities in its gas distribution. Likewise, it is the combination of these two gas phases that will give the 
tightest constraints on the galaxy's history. Additional matched resolution HI and CO maps of Virgo spirals would be 
very helpful for understanding the evolution of the ISM in mildly disturbed galaxies. 

6. SUMMARY 

Images of HI and CO emission at 15" resolution in the Virgo Cluster spiral NGC 4647 show that the neutral gas 
is strongly lopsided. A spiral arm or a segment of a ring on the east side of the galaxy shows up as enhanced HI 
column densities and a plateau in the molecular gas distribution. The west side shows a very steep dropoff in molecular 
content but atomic gas dribbles out to large radii. Thus, at radii larger than 30" (2.5 kpc) the ISM is much more 
highly molecular on the east side of the galaxy than on the west side. As a result, total gas surface densities are twice 
as high on the east side of the galaxy as on the west side. The galaxy also shows a kinematic asymmetry; the rotation 
curve flattens on the west side of the galaxy but continues to rise on the east side. The velocity difference between the 
two sides is nearly 30% at the edge of the HI distribution. 

An analysis of the surface densities of atomic gas, molecular gas, and the stellar disk suggests that the enhanced 
molecular fraction on the east side can be explained if the interstellar pressure is higher there than in the west. The 
implied pressure difference is about a factor of two (2.4x10'* cm~'^ K) in the outer regions of the disk. We conclude 
that the pressure model of Blitz & Rosolowskv (2004) is probably consistent with the gas distribution in NGC 4647 
since there are several different plausible mechanisms that could provide such a pressure difference. 

We discuss possible causes of the kinematic and gas distribution asymmetries in NGC 4647, focusing mainly on 
ram pressure and on asymmetries in the gravitational potential. The facts that (1) the gas distribution asymmetry 
is strong in the molecular gas but weak in the atomic gas and (2) the galaxy kinematics are less disturbed on the 
side where the total gas density is a factor of two lower lead us to suggest that ram pressure by itself may not be a 
satisfactory explanation for the gas pileup on the east side. In addition, attributing the pressure difference directly 
to NGC 4649 requires the two galaxies to be so close that tidal forces should destroy the disk of NGC 4647. Thus 
we favor attributing both the kinematic and the morphological asymmetries to a lopsided gravitational potential that 
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may or may not be related to the presence of NGC 4649. Gas settling into a potential well might explain the larger 
velocities on the east side, and shocks could convert the gas to molecular form. 

We also note several remaining questions that are raised by these results and suggest directions for future work. For 
example, observations of the stellar kinematics of NGC 4649 should help test for a gravitational interaction between the 
two galaxies. Stellar kinematics of NGC 4647 could distinguish whether the forces disturbing its gas are gravitational 
or hydrodynamical in origin. Hydrodynamic simulations could test whether the kinematic asymmetry drives the large- 
scale variation in the molecular fraction of the gas. A better understanding of the cause(s) of the disturbance will help 
to predict how long the disturbances will remain detectable. 

The comparison of the low density atomic medium and the high density molecular medium is vital to obtaining a 
complete picture of the ISM in this mildly disturbed galaxy. Furthermore, all indications are that a large fraction 
of spiral galaxies (both in clusters like Virgo and in the field) have kinematic disturbances like those in NGC 4647. 
Additional comparisons of atomic and molecular gas in mildly disturbed galaxies should offer important insights into 
the evolution of the molecular phase as galaxies settle into groups and clusters. 

Finally, star formation activity (traced by Ha emission) follows the CO intensity very well on scales larger than 
about a kpc. By this measure the local star formation efficiency is constant throughout the galaxy to within factors 
of two to three. 

This paper made use of the Digitized Sky Surveys, produced at the Space Telescope Science Institute under U.S. 
Government grant NAG W-2166. The images of these surveys are based on photographic data obtained using the 
Oschin Schmidt Telescope on Palomar Mountain and the UK Schmidt Telescope. The plates were processed into 
the present compressed digital form with the permission of these institutions. Ha and R images of NGC 4647 were 
obtained through the NASA/IPAC Extragalactic Database (NED), which is operated by the Jet Propulsion Laboratory, 
California Institute of Technology, under contract with the National Aeronautics and Space Administration. This 
work has been supported by National Science Foundation grants AST 0074709 and AST 0502605. Thanks also to M. 
Cappellari for assistance with his MCE decomposition routines. 
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Fig. 1. — Mosaiccd field of view for the CO data. The nominal sensitivity to a point source (tlic primary beam, or tlic gain of the image) 
is 0.99 within tlic inner contour and falls to 0.5 at the outer contour. Contours arc superposed on an image of the system from the blue 
plates of the Digitized Sky Survey. The primary beam of the VLA at 20 cm is much larger; the sensitivity of the HI data falls to 0.5 at a 
distance of 15.5' from the pointing center and it is effectively uniform over the region shown here. 
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Fig. 2. — HI and CO spectra of NGC 4647. Atomic gas is indicated with the solid line and the scale on the left side; molecular gas is 
indicated in the dashed line and the right side scale. Spectra are constructed by first defining the spatial region of emission in an integrated 
intensity image and summing over the same spatial region for all channels. Velocities are heliocentric, in the optical definition. 
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Fig. 3. — HI, CO, stellar, and Ha emission from NGC 4647. Top right: HI column density from the high resolution (15") data 
cube. Contours are 0.5,1,2,3,5,7,9, and 12 xlO'^'' cm~-^. Top left: CO column density. Assuming a H2/CO conversion factor of 2.0x10^" 
cm"'^ (Kkms~^)~^, a CO intensity of 7.2 Jy kms~^ is equivalent to an H2 column density of 6.7x10'^'' cm~'^. Contours are (—2, — 1, 
1, 2, 3, 4, 5, 6, 7) X 6.7 x 10^*^ cm~^. The gain=99% and 50% contours are indicated as well with solid black lines; note that the noise level 
increases dramatically in the regions beyond gain=99% because of large drops in the effective time-on-source at those distant positions. 
Ellip ses in the upper right corners of the frames indicate the spatial resolution. Bottom left: broadband optical image (nearly Mould R) 
from lKoopmann. Kennev. He Younel 12001]) . A smooth elliptical model of NGC 4649 has been subtracted from the im age, and contours are 
spaced by factors of two. Bottom right: continuum-subtracted Ho image from lKoopmann. Kennev. Younel 120011) . 
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Fig. 4. — Atomic gas column densities, averaged in semi-elliptical annuli. The east (receding) half of the galaxy is shown in squares, 
and the west (approaching) half is in triangles. The column densities are deprojected to face-on and are multiplied by a factor of 1.36 to 
account for the presence of helium. The larger error bar in the upper right corner indicates a typical magnitude of the dispersion about the 
mean in each annular semi-ellipse; the smaller error bax indicates a typical magnitude of the uncertainty in the estimated mean for each 
annulus. These error bars are ±1(T, i.e. they are 2a in length. 
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Fig. 5. — Atomic and molecular column densities, averaged in semi-elliptical annuli. The same HI column density data that are shown 
in Figure|^are reproduced here in open symbols, and the molecular column densities are in filled symbols. In both cases, the east half of 
the galaxy is shown with squares and the west half with triangles. Typical 2a error bars for the molecular column densities are indicated; 
see the caption to Figure|l]for explanation. 
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Fig. 6. — HI velocity field as indicated by the fitted center velocities at each position. The wedge on the right hand side of the figure 
indicates the greyscale values, and contours run from 1340 km s~^ to 1520 km s~^ in steps of 20 km s~^. The cross in the center of the 
figure indicates the adopted center of the galaxy (the position of the nucleus in the R-band image) . The ellipse in the bottom right corner 
shows the resolution of these data. 
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Fig. 7. Individual channel maps. HI emission is in greyscale. Contours show CO emission at —4, —2, 2, 4, 6, and 8 times the rms 
noise level (100 mjy b^^). The dark ellipse in every frame indicates the angular resolution of both the HI and CO, and the velocity of each 
channel (km s^^) is indicated in the top left corner. A single large section of a roughly hexagonal contour indicates the points where the 
gain of the CO image drops to 0.99, and small sections of the gain = 0.5 contour are visible in the extreme corners. 
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Fig. 7. — Continued. 



19 




RA offset (arcsec; J2000) 

Fig. 8. — Major axis position-velocity diagram. HI contours are in red, at levels of —3, —2, 2, 3, 5, 7, 9, and 11 times the rms noise 
level (0.35 mjyh~^). CO contours are in blue, at levels of —3, —2, 2,3,5,7, and 9 times the rms noise level (100 mjyb~^). The cross in 
the lower left corner indicates the velocity resolution (channel width) and the spatial resolution (FWHM of the CO beam; the resolution 
of the HI data is very similar). Other crosses indicate the rotation curve derived from fitted tilted rings to the two halves of the galaxy, as 
indicated in the text and as shown again in Figure |3] The noise level of the CO data increases towards the right-hand side of the figure 
because the sensitivity drops off toward the edges of the mosaiced field of view. 
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Fig. 9. — HI rotation curve for the two sides fitted independently. Open circles show the east side of the galaxy and filled circles show 
the west side. 
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Fig. 11. — H2/H I column density ratio vs. pressure. The approaching (west) side of NGC 4647 is shown with filled circles and the 
receding (east) side with open circles. The same relationships for 13 galaxies in Blitz & Rosolowskv (200^ are shown as gray lines. For 
high pressures (near the center of the galaxy), the molecular gas fraction agrees well in the two halves. The relationship appears to be the 
same as for other galaxies. However, at low apparent pressures, the column density ratios diverge in the two halves. 



